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ABSTRACT 

Previous  work  on  diffraction  by  a  wire  gi'ating  is  briefly 
reviewed.  The  problem  of  diffraction  by  a  finite  grating  consisting 
of  perfectly  conducting  wires  of  arbitrary  cross-section  is  foraaa- 
latod  in  terms  of  characteristic  plane  waves  corresponding  to  the 
various  order  spectra  defined  in  optics.   Scattering  matrix  ele- 
ments are  expressed  as  stationary  functionala  of  current  distri- 
bution on  the  grating  wires,  for  the  incident  wave  falling  at  right 
angles  to  the  grating  elements  and  polarized  either  parallel  or  per- 
pendicular to  them.   Similar  expressions  are  obtained  for  the  im- 
pedance matrix  elenents;   thesis  are  evaluated  for  the  thin  wire 
grating,  and  an  equivalent  network  constructed. 
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1.  Introduction 

The  wire  grating  la  coamonly  uaed  to  meaaura  wavelength  of  infrea-od 
radiation.  Vfith  the  Increased  interest  in  the  frequency  range  bridging  the  far 
infrared  and  microwaves,  it  la  important  to  Investigate  the  hehavlor  of  the 
grating  when  the  wavelength  is  of  the  sane  order  of  magnitude  as  the  grating 
constant. 

In  the  infrared  region  the  grating  constant  is  lauch  greater  than  a 
wavelength.  In  this  case  one  can  assume  that  there  is  no  interaction  "between 
the  wires  of  the  grating.  The  diffraction  field  is  taken  as  the  sum  of  diffrac- 
tion fields  of  the  individual  wires  of  the  grating,  the  field  of  any  one  wire 
calculated  neglecting  the  other  wires  .  The  grating  is  assumed  to  consist  of 
a  set  of  infinite,  parallel,  equidistant,  perfectly  condacting  circular  cylinders. 
The  get  may  "bo  taken  as  finite,  shov/ing  the  hackground  in  addition  to  the  various 

order  spectra.  Only  the  infinite  grating  will  he  considered  here.  This  approach 

2 
was  extended  somewhat  by  Twersky  who  included  the  first  order  interaction  effects 

between  neighboring  elements  in  the  case  of  an  acoustic  wave  incident  on  a  grat- 
ing of  perfectly  rigid  circular  cylinders, 

Ignatowski  took  a  completely  different  approach  to  the  problem.  In  a 

number  of  papers-^  he  developed  a  method  for  calculating  the  field  from  the  cur- 

k 
rent  distribution.  He  then  took  up  the  problem  of  the  infinite  grating  ,  and 

showed  that  dae  to  the  periodicity  of  the  boundary  conditions  it  is  natural  to 
perform  a  Fourier  analysis  of  the  field,  thus  splitting  it  into  characteristic 
waves.  These  characteristic  waves  are  precisely  the  various  order  spectra 
treated  In  standard  books  on  optics.  The  current  distribution  on  the  wires  Is 
idoutical,  on  neighboring  elements,  except  for  a  constant  phase  difference.  The 
Intensity  and.  phase  of  each  diffracted  wave  can  then  be  calculated  from  the  current 
distribution  on  a  single  wire  and  the  phase  difference.  It  remains  to  determine 
the  current  distribution.  One  can  obtain  good  approximations  to  the  current  dis- 
tribution very  simply  at  wavelengths  much  smaller  or  much  larger  than  wire  thick- 
xxeas.  These  approximations,  for  the  case  of  the  circular  cylinder  can  be  found 


1.  See  for  exsjaple   Schaefer  &   Reiche,  Sin  Beitrag  zur  Theorie  der  Gltterinversion, 
Ann.  der  Phys.  32,  577,  I91O. 

2.  V.  Twersky,  On  the  Tkaory  of  the  Uon-Specular  Eeflection  of  Sound,  PhD.Dlsser- 
tation.  New  York  University,  19^9. 

3.  W.von  Ignatowsky,  Diffraction  ernd  Reflection  Derived  from  Maxwell's  Equations, 
Annalen  der  Physik  23,  875,  1907; 25,99,1903; 26,1031,1908. 

U.  W.von  Ignatowsky,  On  the  Theory  of  the  Grating,  ibid,  U4,  369,  I9IU. 
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In  Ignatowsky's  paper.  Although  'bottor  approximations  to  the  current  distri- 
"bution  caja.  "be  obtained  from  an  integral  equation,  calculations  'become  comher- 
some  and  lengthy. 

There  is  another  method  of  calculating  diffracted  wave  amplitudes 
which  was  introduced  by  Schwinger  and  his  co-workers  at  the  Radiation  Labora- 
tory for  the  solution  of  waveguide  problems.  The  variational  method  has  been 

5  6 
applied  to  a  variety  of  scattering  problems   .  These,  however,  have  been  of 

two  types:  problems  in  which  the  propagation  constant  has  a  continuous  spectrum-', 
,  and  those  in  vdiich  it  has  a  discrete  spectrum  with  a  single  imaginary  value,  the 
rest  real  .  This  last  restriction  corresponds  to  the  asstiaption  of  such  dimen- 
sions of  the  Trfaveguide  that  only  the  dominant  mode  is  propagated  without  attenu- 
ation.  If  this  assumption  is  dropped,  the  waveguide  problem'  is  verj-  closely 
related  to  the  grating  problem  with  no  restriction  of  wavelength. 

Just  as  in  the  case  of  a  waveguide  we  can  choose  between  the  traveling 
wave  representation  of  the  field,  and  calculate  the  scattering  matrix,  or  the 
standing  wave  representation  -  and  calculate  the  impedance  matrix.  A  thorou^ 

discussion  of  the  impedance  representation  of  the  diffraction  grating,  with  the 

g 
restriction  to  a  single  propagated  mode,  can  be  found  in  a  report  by  N.Marcuvitz, 

Although  the  results  of  the  calculations  are  listed  there  for  a  number  of  wire 

cross  sections.   ao  indication  is  given  of  the  method  used.  According  to 

9 
Dr.  Marcuvitz  the  integral  equation  was  solved  directly.-^ 

The  variational  method  has  been  used  to  calculate  the  surface  impedance 
of  a  wire  grating  •   .   In  this  case  the  calculations  give  the  total  power  trans- 
fer through  the  grating,  rather  than  the  intensities  of  the  various  diffracted 
waves.   The  relation  between  surface  impedance  and  an  impedance  matrix  element 
for  the  infinite  grating  is  essentially  the  same  as  that  between  the  total  and 
differential  scattering  cross  section  of  a  body. 

5, Schwinger  &  Levine,  On  the  Theory  of  Diffraction  by  an  Aperture  in  an  Infinite 
Plane  Screen.  Phys.  Res.  jk,   958,U9US;  75,lU23,19U9.   See  also  W.  Sollfrey, 
The  Variational  Solution  of  Scattering  Problems,  Research  Report  31-11,  Mathe- 
matics Research  Group,  H.Y.U. ,  Contract  No.  AF-19(122)-U2.   i^.  Sollfrey  sug- 
gested the  method  used  here  and  in  Ref.  1 ,   and  contributed  in  both  cases  by 
many  helpful  discussions. 

6.D.S.  Saxon,  Hotes  on  Discontinuities  in  Waveguides,  lectures  by  J. S. Schwinger. 

7, J,  Greig,  B,  Fleishman,  and  J.  Shmoys,  unpublished  work. 

S.H.  llarcuvitz.  Impedance  Description  of  Planar  Gratings,  Rep.R-lHO-l+7,PIB-90. 

9. Contract  Ko,  WLEKG,  V/2S-099-2c-ll+6,  Polytechnic  Research  Institute, Brooklyn, N. Y. , 
By  personal  communication. 
10. W.H, Watson,  Physical  Principles  of  Waveguide  Transmission  and  Antenna  Systems, 

Oxford  University  Press,  19^7,  PP.63-65. 
ll.J.W.Miles,  The  Diffraction  of  Plane  Wave  Through  a  Grating,  Quart.  Appl.Math. , 

7.  ^5.  191+9. 
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In  this  report  tho  Infinlto  plane  grating  con3i3tlng  of  parfactly 
conducting  viros  of  arbitrary  croos  ooctlon  will  "be  considered.  The  Incident 
wave  will  'be  aasumed  to  fall  on  the  grating  In  a  direction  perpendicular  to 
the  wire,  and  to  have  arbitrary  polarization.  For  this  case  the  scattering 
matrix  eleaant3  will  be  expressed  as  stationary  functional? of  current  distri- 
bution on  the  grating  wires.  Siallar  expressions  will  be  obtained  for  loped- 
anca  matrix  elenonts.  The  Impedance  matrix  elements  will  be  evaluated  for  the 
thin  wire  grating,  and  an  equivalent  network  constructed.  Calculations  for  more 
Involved  cases  are  In  progress  and  will  be  presented  In  a  future  report. 

2.  Foi'aulatlon  of  the  Problca 

It  will  be  assumed  that  an  electromagnetic  wave  Is  incident  on  the 

grating  In  a  direction  perpendicular  to  the  grating  elements  (see  Fig.  l).  It 

g 
was  shown  by  Marcuvltz  that,  for  arbitrary  polarization  of  the  Incident  wave, 

the  problem  can  be  split  Into  separate  problems  for  each  component  of  polariza- 
tion. 

1.  Electric  field  vector  of  the  Incident  wave  Is  parallel  to  grating 
elements.  In  this  case  the  electric  field  has  only  one  component  E  =  ^(x,z). 
This  component  satisfies  the  wave  equation  ^    p      ■^   i2n/\)  p     =0,  with  (^ 
vanishing  on  the  perfectly  conducting  surface  of  the  wires,  and  satiafying  the 
Sommerfeld  radiation  condition,  except  in  the  direction  of  the  Incident  wave, 

2.  Magnetic  field  vector  of  the  Incident  wave  is  parallel  to  grating 

elements.   In  this  case  the  magnetic  field  has  only  one  component  H  =  (Z'  (x,z). 

p  /        p      y   ' 
This  component  catisfles  the  wave  equation  V  0  +  (2n/X)     «=  0,  with  the 

normal  derivative  of  Ct  VLjaishing  on  the  surface  of  the  wires,  and   satisfying 

the  radiation  condition  except  in  the  direction  of  the  incident  wave, 

3,  Field  Eepreaentation  in  Traveling  Waves 

In  case  of  either  polarization  the  electromagnetic  field  can  be  repre- 
sented as  a  sum  of  an  Infinite  number  of  characteristic  modes.  Let  the  incident 
wave  fall  on  the  grating  in  a  direction  making  an  angle  0  with  the  z-axls  (see 
Fig.  1).  This  means  that  the  x  dependence  of  the  Incident  wave  is  such  that  the 
phase  difference  between  currents  excited  in  neighboring  grating  elements  is 
2na  sin  O/X.  The  set  of  modes  is  then  constructed  by  requiring  the  sane  phase 
difference  between  grating  elements.   It  turns  out  that,  in  deriving  stationary 
expressions,  it  is  also  necessary  to  consider  the  conjugate  problem.  In  which 
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the  phase  difference  is  -2na  sin  O/X.  Thu.8,  taking  the  time  dependence  to  "be 
exp(-icot),  the  x  dependence  mast  "be  exp  1  i(2n  sin  0/\  +  2nn/a)x  and  we  have 
the  cosiplete  traveling  vave  representation  of  ^ 

r    n  p  8     "     exp  i  [sCk  sin  0  +  2rrn/a)x  +  p  K  z] 
(1) 

k  «  2n/X  K^  =  k*^  -  (k  sin  0  +  2nn/a)^ 

where     n-mode  rnimher,   an  aroitrary  integer  or  ssro 

p  *=     -  1»   depending  on  whether  the  wave  moves   in  the 
positive  or  negative  z  direction 

^     8  =     _  1,   the  positive  value  corresponding  to  the  original 
proolem,    eind  the  negative  to  the  conjugate  one. 

Only  for  a  finite  numher  of  these  modes  K     is  real;      for  other  values  of  n  the 

n 

choice  of  the  sign  of  one  radical  for  K  is  dictated  "by  the  requirement  of  finite 

fields  far  away  from  the  grating.  Thus  for  jkzj»  1,  the  attenuated  modes  csin  he 

neglected  and  the  field  consists  of  propagated  modes  only  (see  7ig.  2),  i.e.  those 

satisfying  the  inequality 

(2)  -  (1  "^  8in  0)  a/X  <  n  <  (1  -  sin  0)  a/X 


The  special  case  in  which  either  bound  is  integral  will  not  he  considered.  In 
general,  if  any  one  mode  0  is  incident  on  the  grating,  then  all  other  modes 


a.p.s 
of  the  set  UJ     ,   ,   will  he  excited. 


h,   T  S  Wave  Prohlea. 

4.1  The  Green's  Ponction 

Let  us  define  a  Green's  function  for  the  first  problem  (  U^  =  E   )         the 
field  at   {x,z)  doe  to  a  set  of  current  element  at   (x'  ■•■  na,   z'),-a/2  <  x'    <  a/2, 
of  unit  aniplitude  and  differing  in  phase  by  ska  sin  0  .     The  fanction  satisfies  the 
inhomogeneous  wave  equation     . 

00 

(3)  i\r^  *  k^)  G  ia;   x,z;   x\z*)   =  '^_^     exp(  inska  sin  O)  S  (x-x'-na)      (z-z') 


n=-oo 


12.   See  ref.   6      p.  2,   Section  I. 
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Ihe  Green's  function  can  te  evaluated  lay  expanding  It  into  the   set  of  nodes 
defined  in  (l).     Let 

CU)  G  «=  ^    exp  [si  (k  sin  0  +  2rm/a)  (z-x')lf     (x;  x'.z')  . 

Then 

and 


(W       G^  »  -^^    (k  sin  0  +  2rrn/a)^e35)  bi[(1c  sin  0  +2mi/a)(x-x' )]   f  (z;x',z«). 


(Ub)     (^^+  k^)G  »^  (5^/;>  a^  -^  K^)exppi(k  sin  0  +  2fm/a)(x-x' )1     f  (x;x«,z')  . 

-00  L  i  '^ 

SuTjstltuting  the  expression  from  (3)   for  the  left-hand  member  of  (Ub),   ther.  multi- 
plying both  sides  by  Q^p  [  -si   (k  sin  G  +  2Xn/a)x|  and  integrating  with  respect  to 
X  from  -a/2  to  a/2 

(5)  0^/o)a^  +  ^)  fj,  Cz^^'.s')  «      ^(z-B')/a  \ 

Then,  imposing  the  Sommerfeld  radiation  condition,  suid  requiring  continuity  at 
z  «  s* 

(6)  ^n"  -^n  ^"^  ^^  ^n'  ""'^'^^  ^ 
Integrating  (5)  with  respect  to  z  from  z'-  to  z'+,  and  making  use  of  (6) 

(7)  A^  =  l/2aZ^ 
and  finally 

(S)  G(8;x,z;z',z«)  =  (i/2a)  ^  K^~^  ezpfsi  (k  sin  G  +  2nn/a)(x-x')  +  iK  I  z-» '  ] 

4,2  Integral  Equation        ' 

The  total  field  is  equal  to  the  sum  of  the  incident  field  and  that 
generated  by  the  surface  currents  on  the  wires  of  the  diffraction  grating.  Since 
the  grating  is  infinite,  the  current  distribution  on  each  wire  is  identical,  thus 
pamltting  us  to  suia  the  contributions  of  the  corresponding  current  filaments  on 
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all  viros,  and  then  to  Integrate  around  a  eizigle  wire,  say  one  located  "between 
I  ■  -a/2  and  x  =  a/2.  The  suiniaation  haa  "been  already  performed  in  the  dorlra- 
tlon  of  the  Green's  function.  Let  the  incident  field  consist  of  a  single  inode 
n:p,8.     Denote  "by  I         „(s,e)  the  current  distribution  ex61ted  "by  this  field  in 

the  central  cylinder  of  the  grating.     Then  we  have 


(9) 


/     "  /   inc  "*"  /  8C  "^   An,p,8^^^  t    jG(s;r;r»)  I  (r')  dr' 


where  r  =  x,z  and  the  integral  is  taken  axound  the  wire.  The  current  distrihu- 
tion  is  obtained  "by  requiring  0  to  vanish  on  the  obstacle.  The  integral  equation 
for  I  is  then, 

(10)  for  r  on  wire  0^   „  „(r)  =-  G(B;r;r«)  I^  ^  „(r')  dr« 

/  "iP»»      ^/  "»P»8 

U,3  The  Scattering  Matrix 

Let  the  incident  field  be 

i 

(11^  ?^inc=  Z^^   ^   \^,J  :.,j.,s 

the  sumaation  being  over  all  propagated  aiodes  (n=-i,  .,,J;  p*  -l,+l).  The  scattered 
field  at  large  distance  from  the  grating  also  consists  of  the  propagated  modes,  with 
the  same  value  of  s. 

(12)  fna'     ^       ^  \^s  ^n... 

'  sc       ,      m»P.8  '  m.p.s 

m=-i   p  '^' 

The  scattering  matrix  is  defined  as  the  operator  which  transforms  the  given  vector 
A  into  B 

or,  in  components  _  '• 

B  =  S  A 

(13)  ^ 

B       =  Z_  S    ,     A 
01.  g.  8        mgs  «-nps  nps 

In  order  to  eveiluate  the   elements  of  S,   take,   as   in  Eq.    (9),   a  vector  A  consisting 

of  a  single  cocrponent  and  of  unit  amplitude,    and  evaluate  the  components   of  B.      For 

that  purpose  let  us   investigate  the  asymptotic  value  of  the  Green's  function  far 
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away  from  the  grating.   It  la  clear  that  for  kjx'  -z )  >>1  all  attentiated  mode  coa- 
tributions  in  (8)  can  "be  neglected  and 


(lUa)  s->  00  &(8;a:,a;x«,*')->  (i/2a)  ^L     K     expfsiCk  sin  0  +  2rTm/a)(x-2')+ik  (z-z'J) 


a-> -00  G(B;x,2;x«,z')->(i/2a)  ■>]"  K"^exp  isKk  eln  0  +2raa/a)(x-x«  )+iK  (z'-z)] 


iH-i 


(lUb) 


c^  „-l 


^'/'^^^>T-.-i.^^'^^i2J».W-''^' 


\ 


The  asymptotic  form  of  the  scattered  field  is  then 


5* 


=  +    )  G(a;r;r')  I  (r«)   dr« 

8C  J  n.,P,8 


(15a)     .^co    ^-*    .(1/2.)    i      \'^„.H.,(')/S^..-l,-i">     I„.,.^r.).r. 


y 


(15.)      .^-«     ?S-».(l/2.)     £^>';7„.-l..(r)  fjJ  „.  .!>•;,  ,„,^. 


n.p.a 


(16) 

Similarly 
(17) 


S„       „    <-n,p,8     =  -.  (l/2aK^)  \f  (r«)  I         (r»)  dr» 

n,q,8  •^'  '        B     1/   n,-q,-d  n.p.ft 


S  C-m.-q-s  =  +(i/2aX   )  10  (r» )  I  (r»)  dr« 

Multiplying  both  sides  of  Eq.    (10)  "by  I  ir)       ^  ,   *         *.  ^  *.u       * 

*^  1     X      /     -     m,-g-a       and  Integrating  around  the  wi 

(IS)  -]  (p  (r)     I  (r)  dr  =         I  (r)  G(s;r:r' )  I  (r')dr»dr 

J  /  n,p,d  m.-q.-d  JJ      m.-q.-fi  n,p,^ 
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Combing  (16)  and  (18)  j  ^, 

^   f^r,     i^'^   I    ^^^9      ir)  I     (r)  dr 

and  from  (17 )  ve  get 

,  (6  (r)  I    (r)dr/V    (r')  I      ([rOdr" 

(20)       a  n,-p,-8  <-a,-q,-s   2a 


-fi'^n.p.i'^^^(-«:''"^^m,-q-i'^'^  ^^'<^^ 


rron  (19)  and  (20)  we  can  derive  a  reciprocity  relation.  It  can  "be  easily  verified 
that,  due  to  the  properties  of  the  Green's  function,  the  right-hand  meahers  of  (19) 
and  (20)  are  equal.  Eonce 

(a)  K  S     .       =   K  S 

B  a.q.e  ^a.P.a      a  a.-P.-a  -^  n,q,8 

The  main  advantage  of  expression  (I9)  is  that  it  is  stationary  with 

respect  to  variation  of  the  current  distri'butiono  I     and  I      ahout  the 
*"  ^  n,p,8      m,-q-8 

distributions  satisfying  the  integral  eqoation  (10).  It  can  "be  readily  shown  that 

substituting  1°  +  f  I'  for  I     ,  and  1°     "^  ^  I '  for  I      ,  then  eipand- 

^^  n,p,s   '^  n      Ji»P»8      o.-q-a    m     m,-q-8' 

ing  8  as  a  power  serves  in  £  and  ^  ,  the  coefficients  of  the  linear  terms, €  ando  , 

vanish  if  I      and  1°      satisfy  the  integral  equation  (10).  By  substituting 
n.P.s      m,-q-8 

approximate  current  distributions  in  (19),  one  can  obtain  very  good  estimates  of 

the  elements  of  the  scattering  matrix. 

5,  T  M  V/ave  Problem 

The  second  problem,  i.e.  one  in  which  the  magnetic  field  of  the  incident 
wave  is  parallel  to  the  grating  elements,  differs  from  the  previous  one  only  in 
the  Green's  function  and  the  boundary  condition  to  be  satisfied.  In  the  T  £  prob- 
len  the  magnetic  field  was  in  the  xz-plane  so  that  the  currents  excited  in  the 
grating  elements  were  parallel  to  the  y  azia.  The  Green's  function  was  the  sum  of 
the  radiation  fields  of  a  set  of  current  elements.  The  radiation  field  of  a  single 
current  element,  a  line  current  parallel  to  the  y  axis  is  clearly  cylindrically 
13.  See  vr.  Sollfrey  loc,  cit.  ref.  5,  page  5 
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eymmetrlc.   In  the  present  case  we  deal  not  with  the  radiation  field  of  a  line 
current,  tnit  with  that  of  a  unifora  distrl'bation,  on  a  line  parallel  to  the  y  ajcls, 
of  short  current  elements  perpendicular  to  the  y  axis,  all  in  the  same  direction 
and  in  phase.  Let  g(r,r')  be  the  radiation  field  of  a  single  current  element  in 
the  T  E  problem;   then  ^^>/  g(r,r' ),  the  derivative  of  g  with  respect  to  the 
primed  coordinates,  in  a  direction  perpendicular  both  to  the  y  axis  and  the 
direction  of  the  current  flow,  is  the  radiation  field  in  the  T  M  problem.  Although 
the  above  discussion  applied  to  the  radiation  field  of  a  single  current  element, 
it  is  immediately  applicable  to  the  field  generated  by  currents  on  a  set  of  lines, 
distance  a  apart,  with  the  appropriate  phase  difference  between  successive  lines. 
Thus  ve  have 


(22) 


r  (a;x,!:;x',z»)  =  r-:^  G(s;x,z;x«  ,z«  ) 


where  G  is  defined  by  Eq.  (S).  The  field  is  then 


(23) 


(b    -^  ,        ^f      ^(p  i^)     ^/'a,(s;r;r.)l    ^''^  ^' 

'    '  inc   '  sc   '  n,p,a    J    v  ^tV,^ 


where  G  is  differentiated  in  a  direction  normal  to  the  wire  surface.  Since  the 
normal  derivative  of  tangentlad  magnetic  field  vanishes  on  the  wires,  we  have  the 
integral  eqxiation 


^<h\ 


(2M-)     r  on  wire 


^         (r)       =-)  G        (8;r;r')   I  (r') 


y 


dr' 


By  considering  the  asymptotic  form  of  the  Green's  function  given  in  (lU),  we 
obtain,    as  before,    the   elements   of  the  scattering  matrix 


(25) 

From  (2U)  we  have 

(2b 


S  ^  =  +(i/2a  K 

m,q,8    <-n,p,s  ' 


'  )  l[-^.?^  (r')l  I  (r')  dr' 


jUv  '^'P'^^    1    m.-q-^^  jj    m,-q-^ 


;r)  G         (s:r;r')  I  (r«)   dr'dr 


-1> 

and, finally 


(27)  K  S     ^       =  (-l/2a) 


The  above  expression  has  the  same  stationary  properties  as  (20), 

6.  Standing  VTave  Representation  and  Impedance  Matrix 

Another  form  of  representation  of  the  scattering  properties  of  a  grating 
is  hy  means  of  the  impedance  matrix  rather  than  the  scattering  matrix.  This  is 
the  standing  \ra.ve   representation,  as  opposed  to  the  traveling  wave  representation  used 
in  the  calculation  of  the  scattering  matrix.  For  simplicity  it  will  he  assiimed  that 
the  grating  is  symmetric  ahout  the  zy  plane.  This  implies  that 

( 2S )  S  =  S 

m,q,B  <- n,p,8     m,_q,8  <— n,-p,a 

The  modes  will  nov  he  defined  hy  the  fields 

(29a)  I   n  a     ^  ^a   8^^^  ®^  ^°  ^^  ^^°  ®  "*"  ^^/^^  ^ 

(25h)  V,  U)=   A„  „  Bin  K  z  +  B    cos  K  z  J       ' 

^  n,8     n,8     n     n,s     n  ^ 

Far  away  from  the  grating  the  field  will  approach  asymptotically  the 

stun  of  propagated  modes,  so  that  the  far  field  is  completely  described  by  the  set 

of  coefficients  A  ,  B   ,  n  =  -i,. ,, ,0,... ,  J.  The  far  fields  on  the  two  sides 
ns   ns         •    F  .   w  «» 

of  the  grating  are,  in  general,  different,  so  that  we  shall  use  superscripts  1  and 
2  to  denote  the  coefficients  z  — ^  -co  and  z  ->  *  oo  respectively. 

At  this  point  we  shall  make  one  more  step  in  the  direction  of  equiva- 
lent transmission  line  formalism.  This  formalism  is  natural  for  problems  involv- 
ing transmission  of  electromagnetic  energy  in  a  finite  set  of  modes.  The  concepts 
of  voltage  and  current  follow  directly  from  transverse  electric  and  magnetic  fields. 
The  transmission  line  equations  are   then  a  consequence  of  Maxwell's  equations 
directly.  This  formalism  applied  to  the  problem  of  the  diffraction  grating  was  dis- 
cussed in  detail  by  Karcuvitz  ,  and  only  the  bare  essentials  will  be  restated  here. 


■«■'•■   I 
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As  Isefore  we  distinguish  tetween  the  two  prohlems  arising  from  differ- 
ent polarizationB   of  the  incident  wave.     In  the  T  E  we  deal  with  an  electronag- 
netlc  field  having  components  3       H   ,    and  H     only.      The  relevant  Maxwell's  eqxia- 
tious,   with  time  dependence  exp   ^icjt }  assiimed,   are 


(30a)  (iic/v^)H^  =  -   |~5L 

(30b)  (iWa)H^   »  1^ 


(30c)        -i^^^j 


3x 

^   -  ill 


Now  S  and  H  ,  the  transverse  fields,  maat  ha/e  the  x   dependence  given  in  (29); 
y     X 

In  the  n   mode 

(31a)  E  =  V    (z)  exp  is  (k  sin  ©  +  2nn/a)  x 

/    n,8  ' 

(31h)  H  =  -I    (z)  exp  is  (k  sin  0  +  2nn/a)  x 

and,  from  (30b) 

(31c)  H  =  >^  8(sin  0  +  2rTn/ka)  exp  is  (k  sin  0  +2nn/a)  i 

Substituting  (31)  In  (3O),  we  obtain  the  transmission  line  equations 

(32a)      fj    V^^^  (X)  =  (ik/r^)  I^^^  (z)  S      (iK^/Y^)  I^^^  (z) 

(320)     -fr  la.s  (^^  =  (yk)  [k^  -  (k  sin  0  +2nn/a)^J V^^^Cz)  =  i^'^J^.^^i^) 

\    2  ?!  ^/2 

(32c)  X^  »  I  k^  -  (k  sin  0  +  2nn/a) 

(32d)       '        ■  ^a  'HV^ 
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The  field  In  a  given  mode  is  then  completely  specified  \ia.en.  1         (z)  and  I    (z) 

n,  8         n,s 

are  given  for  sone  value  of  z  =  z*.   Although  z'  is  arbitrary,  its  choice  "being 

deterained  by  convenience,  let  us  take  here  for  definiteness  z'  =  0  for  the 

representation  of  fields  on  both  sides  of  the  grating.  Let  V  ^  =  V  P  (0)=  I  ^   (0)' 

_^  _  n.s         n,3  n.s 

and  define  vectors     V     =    j  V  ■^    (  and  I     =    |  I  ^   i    ;      the  problem  now  reduces   to 
finding  either  the  matrix  Z  or  its   inverse,  I, 

(33a)  .  V     =     Z     I 

(33^)  f    =     Y     V 

In  special  cases  Z  ©r  Y  may  be  sin^lar  thus  making  the  inverse  representation 
impossible. 

Since  the  grating  is  assumed  to  be  symmetric  about   the  z  =  0  plane, 
the  problem  can  be  simplified.     The  matrix  Y  (or  Z),   whose  rank  is   twice  the 

number  of  propagated  modes,  r,  can  be  decomposed  into  four  matrices  of  rank  V. 
Instead  of  using  single  vectors  V  and  I,  introduce  voltage  and  current  vectors 
on  two   sides  of  the  grating 


i^     =rvM       and  I^     =    (iP  ?. 


then 


(33a)  f-     =     Y^  7'     +    Y^2  ^ 

(33b)  I^     =     Y^^  7'     +     Y^^  V^ 

-'^  S  S  6 

=11    =22 
Due  to  the  symmetry  of  the  grating  Y    =  Y   .  Also,  from  conservation  of  energy 


=12    =21 

Y^   =  Y   ,  We  shall  now  consider  two  cases 


1.  The  field  is  even,  V  =  V^  ,  I '  =1^ 

'88*88 

{■^)  ■/  I®     =   (Y^   +  Y^)     V* 

8  8 


The  problem  is  reduced  to  the  half  space  ziO  and  we  have  the  additional 
boundary  condition    3^    =  0,   z=0. 


d  z 
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2.  Tha  field  is  odd.  V^  -  -  V^  ,1^  =  -  I^ 

8       0      Q       a 


7^ 

=  z^i. 

*Z^ 

'-l^ 

8 

8 

8 

(35)  l^  -  U^  -  l"")  v^ 

In  this  case  the  additional  boundary  condition  is  (P  =  0,  a  =  0 
Siiailarly,  we  hare,  in  terns  of  tha  impedance  matrix 

(33c) 

(33d)  v^  =z^i^  +z^  il 

(3Ua)        ■       8    ^        '  8    . 

(35a)  V°  =(Z^^-Z^)I^ 

We  adso  have  the  charaxsteristic  modes  in  the  standing  wave  representa- 
tion, with  the  Eidditional  hoandary  condition  at  z  =  0. 


In  what  follows  the  tern  "incident  field"  will  he  \inderstood  to  denote  the  imperturhed 

0 

,8   • 


8tandin<5  wave  pattern  given  by  a  linear  ooiahlnat ion  of  0      or  of  0 
7.  Integral  Squation  and  Stationary  rorm  of  Impedance  Matrix  Elements 


We  now  proceed  to  derive  the  Groan's  functions  for  the  even  and  odd  fields, 
and  their  asymntotic  forma.  Since  the  Green's  function  for  free  space  was  derived 
previously,  the  additional  boundary  condition  can  be  satisfied  by  the  image  method, 
For  the  even  case  we  have 

G®C8;r;r«)  =  G(s;x,2;x',z' )  +  G(s;x,«;x« ,-z' )  » 

=  (i/2a)  ^^   K^-""  exp  si  (k  8ln  0  +  2Tm/a)(x-x«)  j  (^  ^\V'^'\   *  exp  iK^Jz+B'l) 


O=-00 
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80,  la  the  region  vmdor  considaration,   z  ^  0 

CO  .  (. 

(39a)  G®(3:r;r«)=(l/a)2I  \     «3q)  8i(k9iix-0+  STa/a)(x-x') 


eip(-iK  z)co6  K  z'    z  5  z* 


(3913)       ■  (l/a)^K~^  exp  [8i(kBlnO+  2rca/a)(x-x» )  J  exp  (-iZ^z»)cos  Zz       z  >  z« 
-co 


The  aayniptotic  form  is  ottainad  from  (39a) 


(UO) 
z-> 


-00  G®(«;r;r»>»(i/2a)^K''^expr8i(koiiiO  *  2na/a)(x-x' )J  fo2p(iSj^a">»'o^(-lE^z')|  oxp(-iZj) 


and  for  the  odd  case 


(Ul)  G°(a;r;r«)  =  G(8;x,z;x,z)  -  G(8;x;z:x,-z)  ■ 

«»(i/2a)^K~"'"exp[si(k  sin  0  +  2nm/a)(x-x« )!(  exp  IZ^^I  z-z«j    -€Kp  iK^  |z+z«l  J 

(U2a)  G°  =  (-l/a)2Z  K~^«3qp[3i(k  ^^^  ®  **■  2nm/a)(x-x' )]    exp(-iK  z)sin  K  z»        z   <  z« 

(U2b)  G*^  =  (-1/a)  ^  K~^e:q)  r8i(k  sin  ©  +  2nm/a)(x-x« )]  exp(-ii:^z«)  sin  K^z     z   >  z« 

and  the  asyiflptotic  form  is 

G°->  (i/2a)^K~^expr8i(k  sin  0  +  2nm/a)Cx-x«  )]/exp  iK^(i«-z)-exp  iK^(-z»-z)J 

(U3)  "  _, 

-  (i/2a)  ^K;^yi^^.,^Jr)[/^^,,^^r.)  -  ^^^^^Jr^)\ 

Even  Case 


In  the  even  case  we  have 


f  '■   <i^M  +  jaU-.T;!')  l^{l>)     dp" 
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Let 


c^)     f.M  -  ^  .,  f:^, 


Porthermore,  daa  to  llaearity  of  the  problem, the  ourreat  I  can  "be  roBolved  Into 

o 


eonrponents  dspaadlag  on  propagated  mode  intensities 

Let  then 

(U6) 

?roQ  (Ui+)  and  (H5),  the  a*^  node  roltage  is 


\M    ■    ^    C    Sn%^) 


{\7)    '  °  (s)    *  a„  MS  i„i  +  (l/aS.)  osjiC-lK  s)    hZ*     °  (r>)  I   (r")  dr> 
(U7a)  7^^^  -  V^;^(o)     -  a^  *  UM,)/)!.!^^. )  I^Cr')  dr* 

rrom  (32a) 


d    „  e 


-1 


1/   (z)     =  (Y  /iK   )  —•  V^*(z)  =   (IT  /K   )  K  a^sin  Z  z  +  a"*^exp(-lK  s)     • 
n,8  n       n     dz     n,d  nana  a  a 


(US) 


r')     I  (r')  dr 


[kS&)  I  ®    =  I   ®  (0)     «  (iY  /aS   ) 

a,  8         a,  8  a'      n 


■] 


(U9) 


The  integral  equation  for  the  pro"blem  is 


Z-n^n!^^'    =-/''°'='^' 


r»)     Ijj(r')  dr' 


where  the  integral  is  taken  along  the  half  wire  lying  in  z  <  0, 
function  can  he  eplit  into  two  parts  as  follows: 


The  Green's 


G®(8:r;r«)=(i/2a)  ji^K^'^'exp   slUsia©  •>  2Tra/a)(x-x« )  [«xpil^  [x-z^i-^  oxp  iK^\z  +  s'  \'\ 

*  1^1  K^^   e^qp  8l(lc  BinO  +2nn/a)(x-x«)\8inK^\z-2«|  ♦  eln  K     \z  ♦  x'\~\  f 

J 
+  (i/2a)^K~^exp  aKkoin  0  +  2nn/a)(x-z> )  [  coa  Kjjj\z-zM   +  cos  Kj^\z+z'  \   J 


(50) 
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where  2L  denotes  suamatlon  over  all  non-propagated  aiodes, 
She  Integral  equation  l^acomes  thon 

^K*  C^/^n>   ffn^'^     ^o^^'^  ^'^  Ky^     =  3/f '(ajnr.)  I^(r.)  dr- 

-i 
or,   from  {hj&) 


(51) 


^    K,B    ^  nM^       =-/r  "   («J-^-')     lo^")  ^'' 


and,  making  use  of  (U6) 

Finally,   since  Y'a  are  artltrary,   one  coefficient  of  each  7  ®    jaast  vanish,   thtla 

th  ^'^ 

giving  us  ain  integral  equation  for  the  n      mode  current  distribution, 

(53)  ^,>)  =-/r 'U;r;-')^n!3("^   ^'      ' 

We  now  obtain  the  adalttance  matrix  elements  by  substituting  , 

(Uto)  into  (USa)  . 

C5U)      .   ,   C  MlVaK^)  ^  V,;,//^_>.)^„>.)dr. 

so  that 

(55)      CV"*V^)„_^^^_,  =(^^Mi//*n>' ',>'»>•>  -^ 

Following?  the  seme  procedure  as  in.  the  derivation  of  stationary  expressions   for 
the  scattering  matrix,  we  multiply  (53)  ^V  ^    ^'^     find  Integrate 
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and  va  obtain 

(I^  +  Y^)       --(in/ak)    ^  ^»-^ </n,+^ J/n,h    i/ n,-h 

a,8;n,s      ^  f C  /)  r\ 

It  follows  that  the  law  of  reciprocity  holds  in  the  usual  circuit  sense 

(5S)  (Y^  ->■   Y^)       =  (Y^  -o-  Y^) 

^  m,8;n,8  n.-s;  m.-s 

Farthermore,  "by  considering  the  unperturbed  field  of  the  form  y>  ^  *  r  *  » 
it  can  "be  shown  that  the  admittance  matriz  elements  are  imaginary.  These  un- 
perturbed fields  are  exactly  the  fields  in  a  rectangular  waveguide  considered  "by 

Odd  CaG3 


0 
Schwingar  . 


In  the  same  manner  ve  evaluate  not  the  admittance  matrix,  but  its 
inverse,  the  impedance  matrix,  for  the  odd  field.  This  is  natural  since,  for 
no  perturbation  grating  vanishing  or  absent,  ve  have 


and     (59) 


Y^.Y^ 


Z^-Z^ 


80  that  in  the  perturbation  theory  used  here  the  inverse  matrices  would  be 
singular  in  this  case* 

VTo  have  then  . 

(60)  ^°{r,  .  2:  a„  )Z^„:,(r) 

-1 

(62)     y°Az)  =  a„  sin  K  z     +     (l/aK   )  exp  (-IK^z)  J   (^    °(r' )     I   (r')  dr' 
(62a)  "  C"Va^°^     =-(VaXn^/^n!-i^'^     I^(r.)dr. 
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(63)  i  ^  r  1 

1°     =  I  °   (0)  «  (Y  /IK  )  4-V  °   (sa         -  -lY Ja  -d/aZ   )  1^    °(;r')  I    (r«)dr' 
Splitting  the  Groon'o  functioa  In  the  saae  way  as   in  the  even  caaa 
G°(8;r:r«)=(i/2a)jX'^n^«^[f^(^   s^  ®  +2nn/a)(x-x' )]  (exp  iK^)z-8'f-exp  iZj  z+iM  ) 

+  i2[^K^^©xprsl(k   sin  9  +2nn/a)(x-x» )  J  (sin  Kjz-z'l-sin  K^|t+z«l    )/ 
(61*) 

(i/2a)X!K~^expr8i(k  sin  0+2nn/a)(x-x« )]    ^cos  K^(z-z«)  -  cos  Z^  (i+s')J 
-t 

The  integral  equation  is  then 

yron  the  alsove  and  (6l),    equating  the  coefficiente   of  I 

(66)  ^^°Jr)  =  -i  Y^/r°(s;r;r')^^^^r')  dr'  ^ 

From  (62a)  j  ^ 

(67)  V^°^  .-(IM^)  Z^  I„_,|  9i,^4r.)^  ^^Jr.)  dr. 

SO  that 

(63)  u'^-.^),_3,„_,=-(lM„)J^44r.);':_^r.)dr. 

But,   from  (66) 
and  finally 
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The  dorlyation  of  elnllar  eqaationa  for  tho  IXipedaiice  caxd  adolttance  matrices  for 
a  T  M  wave  followe  Idontlcally,  and  for  that  roaBon  will  not  bo  presented  here. 

8.  Application  to  tha  Thin  l/iro  Gratipy; 

The  iqpedf.nce  matrix  elements  will  now  he  evaluated  approxinately  for 
the  case  of  a  grating  with  sufficiently  thin  wires.  It  turns  out  that  If  all 
hut  the  most  sigaif leant  term  are  dropped,  it  is  poseihle  to  construct  a  very 
simple  equivalent  circuit.  The  purpose  of  this  section  is  to  ohtain  this  cir- 
cuit rather  than  to  evaluate  the  matrix  elements  as  accurately  as  possihle, 
rurthoraore,  onlj'  the  transverse  electric  field  will  he  considered,  since  the 
diffreiction  effects  in  the  treuasverse  magnetic  case  are  nogligihle. 

In  the  case  of  thin  wires  we  can  use  the  so-called  static  approximation. 
Let  us  first  find  an  asymptotic  form  of  the  Green's  function  0®  for  r  r'.  G®  is 
the  radiation  field  of  a  set  of  line  currents  at  (x*  +  na,z')  and  their  images  in 
the  z=0  plane.  If  we  let  )  r-r'  \  he  small  compared  to  a,  the  contrihutions  of  all 
line  currents,  except  that  at  (x',z')  and  its  image,  can  he  neglected.  Further- 
more, if  k  |r-r'I  «  1,  the  fanction  which  gives  the  radiation  field  of  a  single 
current  filament  can  he  approximated  hy  —  log  kjr-r'J.  Since  P  ®  differs  from 

A 

G  hy  a  hounded  fanction,  this  difference  can  he  neglected  compared  to  the 
logarithmic  term. 


tflro  CroBo-Section 

ns.  3 


Thus  ve  see 
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that  the  integral  Jdr  J^dr"  J^^(r)  P  ®(r,r«)^  ^(r«)  can  ho  replaced 

^^  V    Jc  ^  J^*    3  ^^^   ^°^  lr-r«J  ^^(r«)  ,  ^(r')  defined  on  c  hy  reflection 
1    c^-^c^ 

(see  Pig.  3).  IXie  to  syametry,  we  can  perform  the  reflection  with  respect  to  r  and 
the  integpral  hecomes  "^fj^  ^^)   log|r-r«|  ^^(r')  .  If  we  then  let 

J   ^(r)  ''(/ ■jjd')  =  c[(r)  the  static  charge  distrilmtion  on  a  single  wire  in  infinite 

space,  the  integral  is  proportional  to  the  energy  of  the  charge  dlstrihution.   Using 

W 
the  theorem   that  the  introdaction  of  an  uncharged  condactor  into  an  electro-static 

field  of  fixed  charges  diminishes  the  total  energy  of  the  field,  we  conclude  that 

the  energy  is  hounded  from  ahove  "by  that  of  a  strip  contained  hy  the  wire  and  from 

helow  by  that  of  a  cylinder  completely  enclosing  the  wire.  These  hounds  are  easily 

evaluated;   in  the  limit  of  thin  wire  only  the  logarithmic  term  is  retained,  and 

4  p  2  0  /* 

the  hounds   turn  out  to  he    Q  —  log  kd         and    (^  —  log  kD,  ^  =f  q(r)dr.     Assu-ing 

^1 

that  -T   is  hounded,  the  two  limits  differ  hy  at  most  a  constant,  and  since  we  have 
ct 

heen  neglecting  constants  compared  to  the  logarithm,  take  the  lower  hoxind  as  the 
approximate  value  of  the  integral. 

It  remains  to  evaluate  integrals  of  the  type  j  0  (r)  q(r)dr.  Since 

the  wire  is  thin  9  (r)  on  c,  differs  from  its  value  at  r  ,  the  "center"  of  the 
wire,  hy  a  small  quantity,  and  we  can  take   (r  )  as  a  mean  value.  The  value  of 
the  integral  is  then  Q  ^  (r  ). 
We  have  then 

^'n   12'^^^.^^^   ^  2ak       logkD  ^  2ak  log  kD 

taking  r  =  0  .  , 

lU.  Stratton,  Electromagnetic  Waves,  New  York,  I9I+I,  p.  II7 
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For  thin  wires  the  transverse  electric  fields  in  the  propagated  modes 
aTo  very  nearly  continuous  across  the  S5=0  plane  so  that  (Z   -  Z.  )  =  0,  and 

(72)  I^  +  I^  =  2(Y^  +  Y^)  V 

1    2 
where  the  superscript  of  V  was  dropped  since  V  =  V  .  This  means  that  the  grating 

can  "be  represented  "by  a  shunt  network  (see  Fig.  U)  in  which  the  sh\int  current  is 

given  ty  (in  y)/2ak  log  kD)(X  V^  ).   Such  a  network  is  shown  on  Fig.  3.  The 

transformation  ratios  are  equal,  hut  arhitrary,  and  we  can  set 

(73a)  11=  1   ■ 

The  impedance  is  then  -  '   ,      ; 


\  ; 


■»H*Ji.i.cr«ir 
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Fig.  h 


Equivalent  Hetwork  of  Thin  Wire  Orating 
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